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The effective action for the current and density is shown to satisfy an evolution equation, the 
functional generalization of Callan-Symanzik equation. The solution describes the dependence of the 
one-particle irreducible vertex functions on the strength of the quenched disorder and the annealed 
Coulomb interaction. The result is non-perturbative, no small parameter is assumed. The a.c. 
conductivity is obtained by the numerical solution of the evolution equation on finite lattices in the 
", absence of the Coulomb interaction. The static limit is performed and the conductivity is found to 

■ be vanishing beyond a certain threshold of the impurity strength. 
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I. INTRODUCTION 



The currently used method to obtain observables for disordered systems is to perform a quenched averaging, an 
averaging of the Green functions over random external field [1]. The computational algorithm boils down either to 
the analytical continuation in the number of replicas [2] , or to the introduction of fictious particles related to the real 
ones by super-transformations [3] or to the use of the Keldysh contour in computing loop-integrals [4] . The impact of 
' weak disorder on the conductivity has nicely been captured by resumming the dominant graphs [5] in a self-consistent 
manner [6]. But the whole scheme is based on perturbation expansion and it remains unclear how to treat strongly 
■ disordered systems and understand the transition to localization. This feature appears to be a serious limiting factor 
ON \ in the description of strongly correlated systems, as well, where the interactions are supposed to be non-perturbative. 
The goal of this paper is to demonstrate that methods, well established in other domains of physics, in particular 
the density functional scheme and the functional extension of the renormalization group might be useful to tackle 
this problem. Namely, a new method is proposed for the systemical, non-perturbative computation of the queneched 
' average of the one-particle irreducible (1PI) vertex functions. 

The density functional [7] turned out to be a powerful formalism in handle the ground states. Its virtue is that 
the variational setting gives a large degrees of flexibility for improvements of the density functional. The extension 
to include the current in the variational techniques has already been suggested, as well [8]. This method is improved 
I ■ on the following points. First, the density-current functional is generalized for time-dependent density and current 
configurations in order to encompass the dynamics and to arrive at a closed set of equations to determine this 
functional. Second, a constructive definition of the functional is used which is based on the Legendre transform 
of the generator functional for the connected Green functions for the density and current. In fact, it has been 
realized lately that such an effective action can be identified with the density functional [9,10]. Third, the functional 
generalization [11,12] of the traditional renormalization group procedure [13] is used in the internal space [14] to obtain 
the effective action. The evolution equation describes the changes generated by the gradual increase of the amplitude 
of the fluctuations. The only approximation committed in this scheme is the truncation of the effective action, the 
projection of the functional renormalization group equation into a constrained functional space. One hopes to carry 
out this truncation in a manner which is independent of the amplitude of the fluctuations and the resulting scheme 
will be non-perturbative. We shall apply gradient expansion and retain arbitrarily high orders which seems to be a 
consistent scheme in determing the kinetic transport coefficients. 

The functional renormalization group method has already been applied for disordered systems, namely to the 
pinning of an elastic system in a random potential [15]. A random external potential introduced in an elastic system 
generates infinitely many relevant or marginal coupling constants and their treatment required the functional extension 
of the renormalization group procedure. The U.V. cut-off was gradually decreased and the resulting renormalization 
group equation was solved within the framework of the replica method. In the present work we do not touch the 
cut-off, instead we follow the evolution of the system when a coupling constant is gradually 'turned on'. Such a 
modification of the rules of the renormalization group method allows us to embark directly the issue of quenched 
averaging and offers an alternative to the replica method. 
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The initial condition is usually chosen to be a theory with such parameters which allow an approximative determi- 
nation of the effective action. The evolution equation generates the 'renormalized trajectory' along which the physics 
changes since the quantum/thermal fluctuations are turned on gradually. The integration of the evolution equation 
from the initial condition with small fluctuations to the strongly coupled theory under consideration represents a 
general purpose algorithm to solve strongly interactive systems. 

It is shown below that the quenched average of the generator functional for the connected Green functions for 
the density satisfies a simple evolution equation when the strength of the disorder is gradually changed. This is not 
surprising for weak disorder since the perturbative averaging over the impurity potential introduces correlations which 
are represented by introducing 'interaction vertices' for the density and the corresponding perturbation series involves 
Green functions with increasing number of density insertions. The evolution equation is a differential equation for the 
generator functional for the Green functions which expresses the change induced by the infinitesimal increase of the 
strength of the disorder in terms of functional derivatives with respect to the external source coupled to the density. 
In case of an interactive system the coupling constants of the annealed interactions are evolved, too. The equation is 
obtained without any reference to small parameters in the dynamics and is valid non-perturbatively. The evolution 
equation, being a functional differential equation represents a wonderful mathematical problem. In lacking theorems 
and other general support one has to truncate this equation by projecting it into a restricted functional space. The 
ordinary functions characterizing the functionals in this space, the running coupling strengths in the language of the 
renormalization group, will then be evolved by differential equations. It is rather natural that local functionals are 
easier to handle in a truncation scheme. Therefore the evolution equation is recasted in terms of a local functional, 
the effective action for the density. 

The generator functionals are well suited to the linear response formalism since the introduction of the external 
sources can easily be achieved by functional derivation with respect to appropriately chosen external sources. Therefore 
we include the electric current besides the density in our formalism in order to keep track of issues related to charge 
transport and will use the generator functional for the connected Green functions for the current and density operators 
and their effective action. The formal evolution equation will be presented for the case of electrons with Coulomb 
interaction but the numerical solution will be discussed for the non-interactive case only. 

We start in section II with the introduction of the model and its generator functionals. The average of the generator 
functional for the connected Green functions for the current and the density operators over the impurity field is defined 
in the next step. Its functional Legendre transform, the effective action or the current-density functional is the main 
object of this paper. Connection to the density functional theory is briefly commented, as well. It is furthermore 
mentioned that the generator functional is well suited to the linear response computations except that the causal 
propagators have to be replaced by the retarded or advanced ones in order to implement the physical boundary 
conditions in time. This gives a simple short-cut to re-derive the Kubo expressions for the transport coefficients. 

The evolution equation method is discussed in section III. Though the complete interactive case is covered, the 
equation is studied in a more detailed manner for non-interactive particles only. The functional differential equation 
is approximated by truncating either the powers of the fields or the power of the space-time gradients acting on the 
fields. The former is useful to establish a contact with the resummation methods of the perturbation series, the latter 
is more appropriate for the computation of the transport coefficients. 

Before solving the evolution equation we have to construct the initial condition. This is the subject of section IV. 
The initial condition is defined in the absence of disorder where the effective action is obtained by retaining infinitely 
many orders of the gradient expansion. First we compute the one-loop two-point functions for the current and density 
and construct the generator functional. After that we perform the functional Legendre transformation to find the 
effective action. 

The brief survey of the numerical solution of the evolution equation is presented in section V. The extrapolation 
of the conductivity to the homogeneous case found out to be rather difficult for weak disorder and should be studied 
in a more detailed manner. But strong disorder suppresses the long range modes and the finite size dependence and 
allows to study the conductivity on lattices up to the size 320 3 . The preliminary results presented in this paper show 
a discontinuity developing in the static limit, in particular the d.c. conductivity drops to zero at a certain strength of 
the disorder. It is natural to identify this point with the onset of localization. 

Section VI is for the conclusion. The technical details of the functional calculus, a formal relation between the 
evolution equation and the replica method and the perturbative evaluation of the two-point functions for the current 
and density are presented in the appendices together with the derivation of the Ward identities which are needed to 
assure the finiteness of the d.c. conductivity. 
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II. CURRENT-DENSITY FUNCTIONAL 



The model will be introduced in this section together with the generator functional for its connected Green functions, 
its quenched average and functional Lcgcndrc transform, the effective action. 



A. The model 



We consider spinless electrons propagating on an external electromagnetic potentials V x xt and A x xt and eventually 
interacting with the Coulomb potential. The path integral representation of the vacuum-to-vacuum amplitude is 
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c.f. Appendix A for the notations. We shall use a space-time lattice with spacings a,j — a, j — 1,2,3 and ao = a T . 
There several reasons to introduce lattice regularization. The most important one is that the non-relativistic models 
are usually non-renormalizable and need cut-off for their definition. But even if we are satisfied by considering non- 
interacting electrons, as will be the case in most of this paper, which are U.V. finite the functions appearing in 
the description such as the Green functions develop I.R. singularities. This requires a special care in the numerical 
treatment of the one-loop integral of the evolution equation. The breakdown of this loop integral into a large but 
finite sum can be done in a singularity free manner only if the gauge invariance properties are not spoiled in the finite 
sum. The only known method to achieve this is to use lattice regularization. The regulated action will be 
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in terms of the inverse propagator 
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where jl — a^e^ and the link variables 



U (x) = e -^A^+ev: xt +eu x )^ 



U j (x) = e-* e 



(4) 



(5) 



are used to assure gauge invariance [16]. 

The generator functional for the connected Green functions of local operators O a . x ,Oa, x coupled to the external 
sources <r a . X} ba, x , ol = 0, 1, 2, 3, a = 1, 2, 3, is 
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The operators Oj tX are introduced in order to derive simplify the formal steps in deriving the Ward identities. 
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B. Averaging over the disorder 



Let us suppose that the presence of the impurities can be described by a static potential v x , with Gaussian 
distribution [1] . The quenched average of observables corresponding to a typical static impurity configuration can be 
constructed by means of the generator functional for averaged connected Green functions, 

W[a] = J^Me^^Voh+^^i,--^ (8) 
JV[v)e-^L v * 

Two remarks are in order at this point. First, local impurity potential distribution is assumed for the sake of simplicity 
only, each step below can be reproduced with any other Gaussian distribution. Second, such a simple averaging is 
available for the Green functions of the density only, the operator which couples linearly to the impurity potential. 

Since the impurity potential appears as an additional external source coupled to the density one expects that its 
modification of the dynamics can be obtained by means of the external potential tr . A simple explicit equation which 
realizes this possibility can immediately be obtained by rewriting (8) as 



W[a] = Wo 



5 
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This equation can be used to generate the usual perturbation series [1]. 

Another variant of Eq. (8) which will be useful for the derivation of the evolution equation is established by using 
the 'equation of motion' for the impurity potential. The equation of motion expresses the invariance of the path 
integral with respect to the infinitesimal shift of the integral variable, 
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The finite shift v p — > v p — io u= Q^jT where where T is the Euclidean time extent of the system carried out for the 
Fourier transform gives 

„„ , J CMe -^/ p ("-p-T-o,.=o,-p)K-i CT o,.=o, P ) 

W[a] = * ^ — — (11) 

fV[v]e-^L v * 

where shifted functional is 

Wo[a ,(Ti, ...} = W [a 0iU1 = = Tv,a ,wjto,<ri, ■ ■ ■}■ (12) 

This equation explains the use of the internal space blocking method for the computation of quenched averages. In 
this blocking procedure a parameter which controls the amplitude of the fluctuations, say the mass for relativistic 
scalar theories, is driving the evolution. In the present case <?, the average amplitude square of the impurity field, 
plays a similar role. In fact, the impurity field fluctuations just smear out the generator functional within an interval 
of the size 0{- s /g) in the dependence on the external source. The evolution equation introduced below follows as the 
increase of g gradually 'opens up' a large interval for the source in which the disorder acts. 

The functional W[a] will be useful to express the change of the dynamics under an infinitesimal change of a parameter 
of the dynamics but leads to complications when such an evolution equation is truncated. The complications arise 
from the fact that W[cr] belongs to the class of non-local functionals which are very hard to characterize and handle. 
The way out form this problem is the introduction of the current-density functional, the effective action for the 
composite operators O a as the Legendre transform of W[a], 

SW[a} . 5W[a] 
T[p\ = -W[a\ + a a ■ p a + a & ■ p & , Pa x = - , p & x = — . (13) 

OOu.x OCJa. 



This functional is local and offers useful truncation schemes. 
A relation which will later be used a number of occasions is 

S 2 W[a] S 2 T[p] 



Sx,x"$a,(l- (14) 
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C. Variational method 



The Hohenberg-Kohn theorems [7] assure that (i) the density functional T[p] which gives the ground state energy in 
the sector of the Fock space belonging to a given charge density p x is a well defined functional and (ii) the minimum of 
the density functional is reached at the density profile of the true ground state where its value reproduces the energy 
of the ground state. We generalize the density functional first to include the current T[p] — > r[p M ], where p and pj 
(j = 1,2, 3) is the expectation value of the density and the current, respectively. After that we allow the density and 
the current to become space and time dependent. 

This generalization is realized by the effective action (13). In fact, the inverse transformation is carried out by the 
change of variable 

_ST[pl _ST[p[ 

and this relation indicates that the ground state corresponds to the stationary points of the effective action. The 
transformation (13) can be used to argue that the effective action gives the quantum action (— ih times the logarithm 
of the transition amplitude) for real time or the free energy for imaginary time. Therefore the Euclidean effective 
action is bounded from below and the ground state is found at its minimum. Therefore the effective action plays the 
role of the current-density functional whose minimization leads us to the ground state. 



D. Linear response 



As mentioned in the Introduction, the generator functional is well suited to the linear response formalism. To 
demonstrate this point let us consider the system in the presence of an additional external electromagnetic field 
imposed for t < by making the replacement V t oxt — > V{ — V t cxt + A ,t and A^ xt — > A' t = A| xt + A t in (4) where the 
eventual time dependence arising from the external fields is shown explicitly. For the computation of the conductivity 
3k = &k,lEi, we use the electric field E = — d t A and the interaction Hamiltonian H' = — e J x A k ^ x j k _ x . The ground 
state |0t), following the external fields in an adiabatic manner might show time dependence, as well. The expectation 
value {0t\jk,t,x\0t} of the induced gauge invariant current for t = is 

(0 |jfc,0,x|0 ) = ^<Oo|<o,x+fc C/ fe( a; )V' S ,0,x - ^J,o,x C/ fc*( a 0lk,O,*+k|Oo> 

e 2 

~ e(O o |0 fc ,o,x|Oo) - ^A kfi ^(0 \tpl Ox+k U k (x)ip sfi ^ + V] ;0 ,x C7 fe( a; )V' s ,o,x+fe|Oo) 

- e J y Q{ - y } SAj~y A ^y-~ A ^P ■ ( 16 ) 

Since the bra (0o| is the time reversed version of the ket |0o) in the absence of spin the functional derivative in the 
last equation contains the contributions of the bra and the ket with a relative negative sign which corresponds to the 
commutator between the external source and the observable in the Kubo formula. The time inversion can formally 
be realized in the framework of the path integral formalism by using external field which is symmetrized with respect 
to the time inversion and by taking the complex conjugate of the integrand for t > 0. This latter step brings the 
product of the Heavyside function and the functional derivative into the retarded current-current Green function, 

(0o|j M ,x|0o) = -e 2 [ Dfoi y Ae, y - -A kfi ^~p*, (17) 

J y m 

and we have 

in agreement with the Kubo formula. In the last equation we have the imaginary part of the causal propagator 
because the current operator is hermitian. 
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III. EVOLUTION EQUATION 



We introduce the evolution equation and discuss its structure for non-interactive systems first. The generalization 
which includes interaction, the Coulomb potential in particular is presented briefly later. 

In the traditional renormalization group method one orders the degrees of freedom fa , fa , ■ ■ ■ , 4>n according to their 
momenta, p\ < pi < ■ ■ ■ < Pn and eliminates them one-by one, starting with the perturbative, high momentum 
modes. The corresponding blocking for the action is 



giSjv-i(0i 



-,4>n-i) = J d( j )Ne tS N (4> 1 ,<P2,-,4>N) ^ ( 19 ) 



Such a method may be called external space renormalization group since the momentum is a scale in the external 
space, i.e. in the space-time. As an alternative, one may order the fluctuations according to their amplitude, the scale 
parameter in the internal space of quantum field theory [14]. The evolution from A = to A = oo of the partition 
function 

/A {• X 

#!•••/ ^jve^"^' 02 '-^ (20) 
-A J-X 

takes into account the fluctuations in the order of increasing amplitude. A more practical, smooth cut-off in the 
internal space can be implemented by introducing a suppression term in the action S[<j>] — > S[4>] + Sx[(j>], such as 

Sx[4>]=-\fa4>. (21) 

In relativistic field theory this term generates mass and the realizes the Callan-Symanzik scheme [18]. 
The evolution equation can easiest be derived by starting with the rather trivial identity 

TX Z[J] = T\j *>M«* SW+ *" = \\ l?Me* s M+*"«. (22) 

In the second step the expectation value standing in the right hand side is expressed by acting with the functional 
derivatives with respect to the source j on the generator functional Z[j) and the result is a functional differential 
equation for Z[j]. The harmless looking first step is actually a generalization of the Hellman-Feynman theorem [19] 
for the time evolution operator. Different choices for A yield important, well known equations. The choice of A as a 
field variable 4>{x) leads to the Schwingcr-Dyson equation. When A is chosen to be a parameter of a transformation of 
the field variable which leaves a part or the full action invariant then the resulting equations are the Ward identities. 

The disadvantage of the internal space renormalization group method is that the renormalized trajectory does 
not belong a to fixed physics anymore and the evolution of the control parameter A has little to do with the scale 
dependence in the theory. But the advantage is the flexibility in choosing A. In fact, by identifying A with the Planck 
constant or a gauge coupling one can construct evolution schemes where the semiclassical structure remains invariant 
or the gauge invariance is preserved explicitely, respectively. The solution of the evolution from a perturbative 
initial condition with small amplitude fluctuations to the desired level of fluctuations provides us a systematical 
non-perturbative algorithm to compute irreducible Green functions. 



A. Non-interacting particles 

It was mentioned in Section II B that the effects of the disorder and the external source do are related. Wc follow 
up this remark by showing that the infinitesimal change of the strength of the disorder can be reproduced by acting 
on W[a] with functional derivatives with respect to <r - In order to keep the measure of the strength of the disorder, 
g, explicit in the equations we make the replacement g — > gX and look for the evolution equation in A, 

t[p] = -W[a] 

JPMe-^/p |B ---^°'— |a w?MJ p |,- p -^o, M= o,- p | 2 
2g\ 2 jV[v]e~^ L v ~ 

1 JP[ P ]e^i'p |p -'-^^- 1>l V »[ l T]j2)Me-^J/'J x ^ 
2.9A 2 JV[v]e-*kf- vl jV[v]e-^L v * 
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5 2 W[a] 
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(23) 



where the dot stands for d\. We shall impose the initial conditions at A = 0. It is shown in Appendix B that this 
formalism is consistent with the replica method so long the latter is applicable. 

The functional differential equations are too difficult to handle and the step of central importance is to find a 
truncation scheme which renders this equation useful. We consider two truncations schemes following the Landau- 
Ginzburg strategy of a double expansion in the amplitude and the gradient of the fields. The truncation in the 
amplitude has more pedagogical value, it provides a relation between the resummation of the perturbation series and 
the evolution equation. The truncation in the powers of the momenta, the gradient expansion, seems to be a natural 
approximation scheme for the effective action since the averaging over the impurity field leaves the ground state 
translation invariant and the conductivity, as other transport parameters, are closely related to the Taylor expansion 
of the 1PI Green functions in the external momentum. 



B. Coulomb interaction 



We considered the system of non-interacting electrons so far. But the evolution equation method can be used for 
interactive systems, too. The non-trivial point is that both the annealed and the quenched averages can be computed 
in a parallel manner by generalizing the method of Rcf. [10] for disordered systems. 

The generator functional for the connected Green functions for real time involving the rescaled photon field u — > u/e 

is 

e*Wo[<rj] = J P [ u p^p^t] e i^ t -(G- 1 +A)-^ + i ^2^-9«+iO+^^)-« (24) 

We make the rescaling e 2 — > e 2 A and find the evolution equation for the annealed average 

I 



W [a,j] = -^2^2 / A A u xUy)\x=y 



1 U^-^.A.M], (25) 



2e 2 A 2 V SjSj Sj ' ' Sj 

in the electric charge. The functional we need is effective action corresponding to the quenched average of the 
connected Green functions, 

T[p,w) = -W[a,j}+a a -p a +j-w, Pa ,^, Wx = 5 -^±lA, (26) 

0<7 a ,x 0J X 

where 

y[g|j1= /*-*J.% + ^i Jl (27) 

In order to perform the quenched averaging by evolving the system we make the rescaling g — ► gX in the distribution 
of the impurity potential and find 

t[p,w] = -W[a,j] 

1 JPMe^/ P '"- p ^ g °— l V 1., J ]J p |,- p -l^o,- p p 

v /PM e ~^X>- p -* CT "- pl V >,j] 
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The last line, the piece which generates the Coulomb interaction is represented graphically in Fig. 1. 

The formal difference between the annealed and the quenched averages appears in the manner the disconnected 
parts are treated in the last two terms on the right hand side of the second equation. In order to linearize the right 
hand side in the functionals which facilitates the computation of the quenched average we introduce a new functional 



whose quenched average 
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satisfies the evolution equation 



K\a,3\ = - 



1 fV[v]e-^V v -*-±°— - p|2 ^[a,i]L |»,_ p -* 



T a 0,uj=0,-p\ 



+ 



25A 2 JV[v]e-&L vi 
V J V[v]e'^ l vl K [a + v,^, . . . ,j] 



2A 



S 2 K[a,j] 



2 J z <5cTO,0,z&7'0,0,z 

We write L[p, w] — K[a,j] and the corresponding evolution equation is 
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(29) 



(30) 



(31) 



(32) 
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which is shown graphically in Fig. 2. 
Therefore the (28) can be written as 

g f S 2 W[a] 1 
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(33) 



The system of equations (32) and (33) is closed and describes the evolution of the interactive system. 
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C. 0(p 4 ) truncation 



In the rest of this Section we return to the non-interacting case and consider two different approximations for the 
evolution equation. First we use the ansatz 



T[p] = T + Tap & + ^r a g p a pg + ^r a - b & PaP~ b Pt + ^afiAdP^P'bPs 



(34) 



where the combined 'super-index' a = (a,x a ) has been introduced and summation/integration is assumed for double 
indices. The inverse of 



r ->] 

a,6 L ' ' 



(35) 



where 



(f- 1 ). K = r„ t 

\ /a,o a,b 

on the right hand side can be obtained by the Neumann-series, 



(36) 



r[p] 



f (f-f-(tW [p] + T^[ P ])-F + F- (r< 3 ) [p] + r< 4 > [ P ] ) ■ f ■ (r< 3 ) [ P ] + [ P ] ) • f - ■ ■ ■) 

Jx V ' (°>°> : 



:),(0,0,x) 



(37) 



The identification of the coefficients of the same p powers gives the equations 

r 



i,x),(0,0,x) 



T a,bP&Pb 
^a,b,cPdPbP& 

^a,b,cJP*PbP£Pd 



(0,0,x),(0,0,x) 



^(0,0,x),(0,0,x) 

i/> r(3) M ( o,o 

I |jF.(r(^]-r( 3 )M.f.r< 3 )[p]).F 

F • (2T^ [p]-F- [p] - V™ [p]-F- [p] ■ F ■ [p]) • F 

J f ■ (r^ [ P ]-F- [p] - [p]-F- [ P ] ■ f ■ r< 3 > [ P ] 
-2v^ [ P ]-F- r< 3 > [p] • f • r^ 3 ) [p] + r< 3 ) [p] • f ■ [ p ]-f- [ p ]-f- [ p ])-f 



(0,0,x),(0,0,x) 



(38) 



(0,0,x),(0,0,x) 



shown on Fig. 3. It is easy to understand the structure of these equations. They express the principle that an 
infinitesimal change of the strength of the disorder, g — > g + Sg, can be realized in the order O(Sg) by visiting each 
correlation line of the graphs and change its coefficient, represented by the cross on the correlation lines of the graphs. 
The equations express this process in a hierarchical manner, in terms of the one particle irreducible vertices. 

The resummation of the perturbation scries, implicit in the solution of the evolution equation, can easily be seen 
by iteration, 



(«)[ 
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&7Oa,=O,- p 50O, w =O,-p 



(39) 



where each step increases the order of the interactions/correlations present by one and the integration over A' inserts 
the necessary symmetry factors. In the numerical integration of the evolution equation each step A — > X + SX increases 
the order of the interactions/correlations by one and all graphs are resummed as 5 X — > 0. Naturally this resummation 
is performed in the given truncation of the irreducible vertex functions in the momentum space. In particular, the 
the self-energy insertions, the ladder and crossed diagrams of a particle-hole loop [17,6,21] are resummed among other 
contributions by the first term on the right hand side of the third equation in (38). 
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D. Gradient expansion 



We consider now the ansatz 

r [p] = / \Pa,xV a ,(){Px,-idx)Pl3,x ~ U(p x )\ 

J X 

for the non-interactive case where the space-time derivatives act on the function pp. x only, 



(40) 



(41) 



The combined index A = (vo, v\, 1/2, ^3) defines the partial derivative d A — <9q d" 1 &^ 2 d$ 3 . Since static disorder 
appears through the po dependence in the evolution equation we retain the 0((pj) 2 ) terms only which are needed for 
the conductivity. Therefore we shall use the expression 



m = 



with 



,Po,zT u (po,z, sign(Imid )id , -A)p 0) * 
+ ^po,z^ ts (po,z,sign(Imid )id , -A)d k p k , z + ^Pk,z^ ts (Po,z, sign(Imid )id , -A)d k p , z 
+\p3,z (F T {po,z,sign(Imid )id , -A)Tj. k (d) + T L (p , z ,sign(Imido)id 07 -A)L : j tk (d)) p kiZ - U(p z 



(42) 



U(p) = --T t pj + V(p ) 



(43) 



and 



L(d) = ^, T(d) = l-L(d). 



(44) 



The coefficient functions satisfy the relations T ts (d) = -T ts (-d), T tt (d) = T M (-d), T T (d) = T T (-d), T L (d) = 
T L (—d). The effective action contains derivatives in arbitrary high order the truncation is that they act on a single p 
only. The computation of the second functional derivative used in the evolution equation is sketched in Appendix D. 

The evolution of the functions parameterizing the effective action can the easiest be obtained by writing p a ^ x = 
Pa + &Pa,x with J 5p a ,x = and expanding T[p + 5p] in Sp up to terms G(Sp 2 ). We shall use the notation = 
pod^ + kd^ with the property d Po V{ = Vj+1, the 7-functions 7*' = r* + V?r"(p), 7 * s = V^r* s (p), 7^ = T s +T L (p), 
jp = V s + T T (p) , and write the left hand side of the evolution equation as 

r[p] + J |^Po,- P 7 tt ^o,p + i5po,- P j ts (p)pjSp^ p + IjSpj.-p [i T {p)T p + i L {p)L p ] <5p fe , p J (45) 
The right hand side gives 



Jp ^ Jp Jr 



r. I ^po,()(r^ ^ )— 'p,p / $P~i,—r$Pa,rF—p,p+r,—r 0-^ ^ )— p— r,p+r^ —p — r,p,r 0-^ ^ )~P,P ^(^P ); 



(46) 



where fj°' K , T 



-t,t,— r,r 



and T 



are recorded in Appendix D. The second term resums the impurity correlations 



attached to the particle-hole loop, the self-energy insertions together with the vertex corrections. The third term 
contributes at vanishing energy only and will be neglected. 
The 0(6 p°) part of the evolution equation is 



fj= £ [ r-i(o,oh = £ [ 7 p 

2i P )P 2i p 7«7 P L +P 2 (7^) 



(47) 
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The 0(Sp 2 ) part gives 



ir 


9 ~f 


K 


7? + a p 2 o7 f](7 P L ) 2 


2 J P 


[7^7 P L + P 2 (7* s ) 2 ] 2 


it s = 

I r 


9 -d 2 

4 po 




f (7 P L ) 2 

/ P [7^7 P L + P 2 (7^) 2 ] 2 


iT = 




£ 1 


(7 P L ) 2 




3 [ 7 «7p L +P 2 (7^) 2 ] 2 


i L r = 




(7 P L ) 2 





(48) 



One can verify as a consistency check that the O(r ) parts of the equations for j L and -f T agree with (47). Eqs. (48) 
have the form of a one-dimensional diffusion equation with (po, g) as 'space-time' coordinates with a 'space'-dependent 
diffusion constant given by the one-loop integrals. 

The evolution equations need two boundary conditions in p . One can be obtained by noting that the Green 
functions D^ ,v are symmetrical under charge conjugation p — ► 6h 2 /ma 2 — p. Therefore d Po "f = will be chosen as 
a boundary condition at half-filling, p = 3h 2 /ma 2 , po — l/2a 3 and the density will be restricted into the interval 
< Po < l/2a 3 . The other boundary condition is imposed at the lowest density point which is chosen at sufficiently 
small chemical potential in such a manner that only one particle is found in the system. By assuming that the Green 
functions increase linearly with the density in this low density regime the boundary condition d 2 a j = 27/pg will be 
imposed. 

IV. INITIAL CONDITIONS 

After having obtained the evolution equations (47) we need their initial conditions, imposed in the perturbative 
regime. One might choose either A = and use the free effective action or take a small but non- vanishing A and 
compute the initial conditions perturbatively. For a computation carried out at vanishing frequency and on infinitely 
large system one has to start with small, non-vanishing strength of disorder. The reason is that we integrate a 
differential equation in this scheme and have to assume that the A-dependence is differentiable. But non-interacting 
particles move ballistically in the absence of disorder and an arbitrarily weak disorder already induces important 
changes in the transport coefficients. In fact, the system becomes localised in one or two dimensions at arbitrarily 
weak disorder. Sufficiently dilute three dimensional systems should be localised at weak impurity field and the 
construction of the initial conditions at weak disorder is rather non-trivial [20] . 

In a more formal language one may say that the difference between the initial conditions imposed at A = and A w 
lies in the way time reversal symmetry is broken. In order to arrive at finite conductivity one needs the breakdown 
of the time reversal and the space-translation invariance. Neither of these breakdowns takes place in the absence of 
interactions between the electrons and their environment. The interaction with static impurities breaks translation 
invariance. The interactions break the time reversal invariance dynamically in infinite volume. The advantage of 
the perturbative initial conditions imposed at small but non- vanishing A is that they contain the desired symmetry 
breaking pattern. 

But we have to bear in mind that the numerical solution of the evolution equation can be obtained for a finite 
system only. Furthermore, in order to avoid the singular chemical potential dependence in finite systems it is necessary 
to perform the computation at finite frequency and to study the limit u> — > numerically. The A-dependence of the 
Green functions of a finite system or at finite frequency is regular at A = and no difficulty is expected during the 
integration of the evolution equation. Therefore the initial conditions will be constructed in the absence of impurities, 
at A = 0. The time reversal invariance is then broken externally, by the non- vanishing frequency and the transition 
from the ballistic to the diffusive regime will be induced smoothly during the evolution. 

It is worthwhile noting that the issue of localization of non-interacting particles is far from being a one-particle 
problem and the density-current functional is a highly non-trivial object even for non-interacting electrons without 
disorder, at A = 0. This is because the Pauli blocking of the occupied states represents a strong topological many- 
particle interaction which obeys no small parameter. This circumstance is crucial in expecting a singularity in the 
ground state in function of the disorder strength. 
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A. Current and density Green functions 



The first step is the construction of the generator functional for the connected Green functions of the density and 
the current. For this end we write a a ^ x — a a + 5a a , x where J 5a a ^ x — and find the following generator functional 

W[a + 5a] = -iTiTrlog [G" 1 + <fr + 5<f\ = ftTrlogL*- 1 + p* ■ 5a - ^Sa ■ D ■ 5a + 0(5a 3 ) (49) 

where D^ 1 = G^ 1 +<fr, p* az = —ihTi [D ■ O a . z ] = p* , and D( a . x ).(p, y ) — —ihTr [D ■ O a . x ■ D ■ Op iX \. Due to rotational 
invariance one has the form 



L>"(w,q 2 ) iqe-i a iD ts (w,q 2 ; 



jqe* aq D ts (w,q 2 



i> fc (w,q) 



(50) 



with D^ k (u,q) = T^ fc (q)i) T (w,q 2 ) + J L^' c (q)i) L (w,q 2 ), q= fsin^, L = q®q/q 2 and T =1 — L. But on a lattice 
with finite lattice spacing there is no rotational symmetry and the transverse piece becomes more involved. This 
complication will be neglected below. A more detailed expression for the propagator is presented in Appendix 
E. 



B. Variable transformation 



The second step in constructing the density-current functional for the initial condition is the change of variables 
a — > p. For this end we have to invert the relation 



Pa.x = (p a + 5p a . x ) = -ihTv 



D- 



^ OL.X 



= ihTi [D ■ O a , x ] + zftTr [D ■ 5<fr ■ D ■ O a , x ] + 0(p 2 ). (51) 



We write p + po — > p set pj = a 3 ; = in the absence of Coulomb interaction and find 

5a = -I)- 1 ■ 5p 



(52) 



C. Effective action 



Finally, the effective action can be written as 



T[p] = ihTilogD- 1 + a ■ p-^5p- D- 1 ■ 5p + 0(5p 3 ). (53) 



The identification of the part O(5p ) yields 



ifi 

U =-a aPa -yTr\ogD- 1 . (54) 



The 0(5p 2 ) order, the second functional derivative gives 



i^ai^ts Ti t + Li l) \ i^ei^D*' TD T + LD L J ' \ _ zqe ^q^ J_ r + ^-L 



(55) 



where 

#(c,q 2 ) = D«(o; ) q 2 )D i (a; ) q 2 ) + q 2 p ts ( W , q 2 )) 2 . (56) 

Therefore 
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y V q 2 ) = - ^ ( "' q2) 

7>,q 2 )' 



D T (co,q 2 ) 



Li 2, D u (oj, q 2 ) 
7 (w,q ) - 



£(^q 2 ) 

*>,q 2 ) 
D(u;,q 2 ) ' 



7 4 V) = ^^. (57) 



The Euclidean Ward identities require 



w \ m 

and 



D ts (ui, q) = -- (^ fe (^,q) + lj**p* I (58) 



D«( w ,q 2 ) = S!^ (W)q 2 )) (59) 



c.f. Appendix C. It is easy to see that these equations restrict 7^ to be momentum independent. 



V. NUMERICAL RESULTS 



There are two points where the numerical integration can be made faster. First, the one-loop integral in the 
evolution equation requires the knowledge of the 7 functions in the whole Brillouin zone. It was assumed that the 7 
functions depend on p 2 , momentum square only. This assumption was found to be acceptable in the initial condition. 
Note in this respect that D L and -y L remain well defined in a non-isotrop Brillouine zone. The other improvement 
on the speed of the algorithm is the runtime adjustment of the step size Ag during the numerical integration. This 
was carried out by keeping the largest relative increment in the differential equations under 1%. Lattices up to size 
320 3 were used but the results presented below correspond to the volume 80 3 unless it is stated otherwise. The units 
a = m = h = 1 were used to express dimensional quantities. The integration of the evolution equation started at 
g = with Ag = 10~ 6 and ended at g ~ 10 5 , where Ag 10 4 was used. 60 division points were used to monitor the 
density dependence of the 7 functions in the interval 2 x 1CT 6 < po < 0.5. The Euclidean frequency was restricted by 
the bounds 10" 5 < u < 335. 

The Ward identities allow us to express the longitudinal conductivity as 

a -^p- = p 2 ) = -™ " 2 L( v (60) 

ie 2 v ' p 2 7"(w,p 2 ) + LU 2 ~f L (uj) 

indicating that the d.c. conductivity in the homogeneous p = limit corresponds to an I.R. cut-off, a 'mass gap' 
in the propagator D ts (uj, p 2 ) and it can only be non- vanishing due to an C(l/p 2 ) singularity in 7" functions 
remains finite as u> — > 0). A rough, qualitative picture of what happens as the disorder is increased is given by the 
analogy with the diffusion equation. The 7- functions of the initial condition have a singular peak at p w 0. Therefore 
the diffusive evolution smears this peak out and tries to 'fill up' the 7 functions at large densities where they take 
smaller values. Since the propagator D u is 0(1/V) for a system of finite volume V the diffusion is rather slow at low 
densities. Therefore the main ^-dependence will be an increase of the 7 functions for intermediate and large densities. 
As this increase diffuses 'up' in the density it increases the 1PI functions and suppresses the conductivity. 

Let us start with the (/-dependence of ct _1 (cj, p 2 ), shown in Fig. 4 (a) at u> = 2 x 10~ 5 and p = Pi = (0, 0, 7r/80) for 
densities 7.6 x 10~ 3 < po < 0.5. The inverse conductivity displays a plateau at weak disorder and starts to increase at 
g w 0.1 only. In order to see clearer what happens at this crossover point Aa _1 (g) = er -1 (g) — cr _1 (0), the contribution 
due to the disorder, is shown in Fig. 4 (b). For some densities one finds a slight drop in er -1 at the beginning of the 
evolution, the figure shows the values with Acr~ 1 (g) > only. We believe that this weak non-monotonic (/-dependence 
is a finite size effect. It is induced by the evolution equation when the density dependence in the initial conditions is 
distorted by a finite size effect. The latter is due to the constraint 0(— Eq) in the integrand of the zero temperature 
Green functions which induces an artificial step function-like chemical potential dependence in the initial conditions, 
(E8). When this finite size effect is ignored the lesson of Fig. 4 (b) is that the change of the nature of a^ 1 at g s=a 1 in 
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Fig. 4 (a) is not a singularity. What happens is that the increase of tr _1 due to disorder becomes comparable with the 
initial value and disorder starts to play a more important role. Act -1 (<?) plotted for different values of the momentum 
p shows similar behavior. One finds the scaling Aa^ 1 (g) g k with k = 1 for g << 1 and k = 1/2 as g >> 1 without 
any non-analicity in the g dependence. The behavior Aa^ 1 (g) ss g characterizes weak disorder. It is easy to see that 
the condition of the simplification that the g-dependence generated by the evolution equations (48) factorizes as an 
overall multiplicative factor is just the asymptotic scaling Act -1 (g) s=s ^fg. The possibility of any sudden transition at 
p = developing in the thermodynamical and the u> — > limit between the weak and strong disorder region is still 
hidden in the momentum dependence, the issue we turn to now. 

In order to explore the I.R. regime the extrapolation of the conductivity to p = was performed by fitting the p 2 
dependence of cr _1 (w,p 2 ) by a polynomial up to p 12 . The result is dominated by finite size effects at weak disorder. 
In fact, g = 0, a perfectly ordered system should have vanishing a^ 1 at p = 0, i.e. long range Green function D ts . In 
more physical terms, the density of state is a sum of Dirac-delta peaks for finite systems and one needs rather large 
volume in order to approach a continuous density of state which is needed for the Bragg reflections on the lattice 
and for finite conductivity. The dependence of the inverse conductivity on the strength of disorder, depicted in Fig. 
5 (a) supports this expectation. Furthermore one can see that at low densities where the density of states is lower 
the conductivity is strongly suppressed and a- 1 is large and positive for any g. When the negative parts of a 
are discarded as in Fig. 5 (b) one sees the sudden drop of the conductivity as the strength of disorder is increased. 
(T _1 is plotted against the density in Fig. 6 showing that the density dependence of the a.c. conductivity is strongly 
suppressed for strong disorder. 

For the identification of the sudden increase of cr _1 with the onset of the localized phase we have to perform the 
limit u) — > 0. The frequency dependence of a^ 1 is shown in Fig. 7 (a). One recovers the usual a w 1/w behavior at 
any momentum p at high frequencies except a short flattening around u> ~ 10 where the disorder decouples. In fact, 
the evolution is suppressed at high frequencies clue to the smallness of propagators in the initial conditions, c.f. (E8), 
and the scaling a w 1/w is recovered at higher frequencies with a proportionality constant given approximately by the 
initial conditions. As we turn towards low frequencies we find two qualitatively different behaviors for the conductivity 
extrapolated to p = 0. At weaker disorder cr -1 decreases and becomes negative, more precisely is dominated by finite 
size effects and should not be taken seriously when computed at the present size lattices. On the contrary, for stronger 
disorder the computation is free of finite size effects and a decreases as u — > 0. Motivated by the Mott form 



with (k, 17) = (2,4), the values of k 77, loq and C were fitted in Eq. (61) for strong disorder. Mott's result alone 
compares rather poorly with the numerical results as shown in Fig. 7 (b). The simplest choice, (k, rj) = (1,0) 
represents a better approximation but the fit with the results (k, 77) = (1.57, 3.18) shows the presence of logarithmic 
corrections and gives a further improvement. 

There is a discontinuity at g = g cr separating finite and diverging values of a~ x at weak and strong disorder, 
respectively. The critical strength of the disorder, g cr can be estimated by finding where tr -1 changes sign first as g 
is decreased, c.f. Fig. 8. This is how the localization-delocalization transition appears in this scheme. 



A new, non-perturbative scheme is presented in this paper for the description of disordered system. The basic idea 
is to express the effects of an infinitesimal increase of the disorder by changing the density. No small parameter is 
required to obtain the evolution equation but truncation is needed to convert it into a useful numerical algorithm. 
We believe that the gradient expansion is a well suited and consistent truncation scheme for the computation of the 
transport coefficients. 

There is a rather detailed space-time picture for weak localization drawn from the partial resummation of the 
perturbation expansion [21]. The mechanism proposed in this paper is more formal, it is a diffusion process where 
the density and disorder strength are considered as space and time variables. The 1PI functions take large values 
at low density in the initial conditions. The conductivity decreases when the disorder is becoming stronger due to 
the diffusion of the low density peak of the 1PI functions towards higher densities. Localization appears as a special, 
correlated increase of the 1PI functions at different momentum values. 

The numerical results obtained in solving the evolution equation indicates that weak localization is rather difficult 
to establish due to strong finite size effects. In fact, the finite conductivity arises from an C(l/p 2 ) I.R. singularity 
of a 1PI function whose reproduction requires unusually large lattices. But strong disorder suppresses the long range 
correlations and renders the numerical solution more reliable. 




(61) 



VI. CONCLUSIONS 



14 



It was found that the conductivity is vanishing as uj — > when the strength of disorder exceeds a certain threshold. 
But the fit of the numerical results yields fractional powers of uj and logw, indicating that several terms contribute 
in the frequency interval considered. It is natural to identify this threshold with the localization transition. 

The results presented here are preliminary, obviously more careful study of the finite size effects is needed to 
establish contact with perturbation expansion and to identify the weak localization regime. Furthermore, a more 
extended frequency interval must be studied in order to identify the low frequency scaling form of the conductivity 
in the localized regime. Finally, the most important step is the inclusion of the Coulomb interaction in the numerical 
solution. 
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(a) 




FIG. 1. Graphs contributing to the evolution equations of the generator functional W[<t,j] in e 2 . The graph (a) and (b) 
correspond to the one-particle irreducible and reducible contributions, respectively. The solid and dashed line represent the 
current-density and photon insertions, respectively and and the circle with n legs denotes the n-point connected Green function. 
The cut stands for the inverse photon propagator. 




(b) 

FIG. 2. Graphs contributions to the evolution equations of L[p, w] in A. The solid and dashed lines represent the propaga- 
tors (8 2 T / SpSp)^ 1 and (S 2 r/8w8w)^ 1 , respectively and the triangular the vertex S 3 T /SpSpSp. The dot is for the functional 
derivatives 5L/8p and 8 2 L/8p8p. The cut stands for the correlation vertex g. 
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(b) 

FIG. 4. a: Inverse conductivity a 1 (g) as the function of the coupling constant g at uj = 2 x 10 5 . Different lines correspond 
to different densities in the interval 7.6 x 1CT 3 < po < 0.5. b: The same as (a) except for Act -1 = a^ 1 (g) — cr~ 1 (0). Both 
a -1 and Act -1 are monotonically decreasing function of the density for g > 1. The negative Act -1 values were ignored in the 
logarithmic plot. 
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(b) 

FIG. 5. (a): Inverse conductivity extrapolated to vanishing momentum at u> = 0.0002 as the function of g for 
0.0171 < po < 0.19. The solid and dotted lines correspond to lattices 80 3 and 160 3 , respectively. The finite size depen- 
dence is stronger for a^ 1 < 0. (b): The extrapolated a' 1 as the function of g for densities 0.0171 < po < 0.5. Only the values 
a -1 > are shown and cr _1 is a monotonically decreasing function of po in this regime. 
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FIG. 6. The positive values of the extrapolated inverse conductivity at a> = 0.00002 as the function of the density for different 
strength of disorder, 180 < g < 4 x 10 4 . a -1 is increasing monotonically with g at high densities. 
















(2,4) 






(1.0) 

















O.OOOl 0.001 0.01 



(b) 



20 



FIG. 7. (a):The solid lines show the positive part of the inverse conductivity extrapolated to vanishing momentum as the 
function of the frequency for po = 0.11 and g — 10, 10 2 , 10 3 , 10 4 , and 10 5 . Lines with symbols correspond to the inverse 
conductivity with momentum p = pi = (0, 0, 7r/80). is a monotonically increasing function of g. (b): The ratio aj 1 /cr^ 1 
plotted as the function of uj in the localised regime, at g — 10 2 , for different choices of (k,t]) shown beside the curves. Notice 
that <r _1 changes by three order of magnitude in the frequency regime shown in this plot. 



100 




0.01 0.1 

FIG. 8. The average strength of the disorder potential in units of the band width, ^/g cr / Ef on 80 3 lattice as the function 
of the density for u> = 0.00512 (heavy line) and 0.00002 (thin line). The dots represent the same quantity obtained on 320 3 
lattice. The finite size dependence is invisible on the lower frequency curve. 



APPENDIX A: NOTATIONS 



We summarize here briefly the notations used in the paper. The space-time integrals and the Fourier transform are 
' l2 it=[, [d d x /'• //,//, /•//. (*' l> ' 

-T/2 



(Al) 

with x = (x, t), p = (p, —it)), and xp = xp — tu). The Dirac-deltas 5 XiV — a~ d a T 8^ v and S p . q — VT8 pq are expressed 
in terms of the Kronccker-deltas Sif,. and 5^„. Identities for the functional derivatives are 

& _ d \ " -i px & _ 1 p-i-P* - = f p-iP x — 

Sf x da*a T f x 2f df p VT^ Sf p J p 5f p ' 

w = VT wr = E * ipx 4r = I ^ vx 4r- ( A2 ) 

ofp d fp x d fx Jx of x 

Greek and Latin indices will be used for the range 0, 1, 2, 3 and 1, 2, 3, respectively. 



APPENDIX B: ANNEALED V.S. QUENCHED EVOLUTION 



It is shown that the annealed and the quenched averages correspond to simple evolution equation. Let us take a 
function f(w) and introduce the generator functional 
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f(W[a]) = 



J V[v}e- s ^f(W [<T-iv}) 
J V[v]e- S -*M 



where 



Sd[v] = -v-K-v. 

We make the replacement v — > v — a (M), v ^> v — ia (E), i.e. use the equation of motion for v to write 

/ V[v]e- S ^ v - 1 ^ f(W a [-iv]) 



f(W[*}) = 



fD[v]e- s *W 



(Bl) 



(B2) 



(B3) 



We therefore have the relations 



2/©[ v ; 



e 2 



(5(7,, 



SW M [ 
5a x 



fT>[v]e- s *W 

JV[v] {(K ■ (v - ia)) x {K ■ (v - ia)) v - K XiV ) e~\^ ) f{w E { _ iy]) 



(B4) 



JV[v]e- 



■iv-K-v 



which allow us to write the evolution equation as 

1 

W[a] = --TrK- 1 ■ 



S 2 W[a] 8W[a] f"(W[a}) SW[a] 



(5(7(5(7 



5a f'(W[a}) 8a 



The simplest equation arises from the assumption 



f"(W[a}) = -f{W[a]) 



(B5) 



(B6) 



where n is a constant. The solution, 

/(*) = {*-/» njS ( B? ) 
corresponds to the quenched average for n = and to n annealed replicas when n^O. The evolution equation reads 



as 



One finds 



2 \ (5(7(5(7 h da da J 



t[a] = -TrK- 1 



SpSp 



(B8) 



(B9) 



in terms of the effective action r[er]. The condition for the applicability of the replica method is the assumption that 
the solution of this differential equation not only converges as n — > but the n-dependence is continuous, i.e. the 
limit agrees with the solution obtained with n = 0. According to Eq. (B8) n should approaches zero at least with 
the inverse volume in the thermodynamical limit. This requirement, well known from elementary considerations, may 
pose problem in phase transitions or when some symmetry is broken spontaneously. 



APPENDIX C: WARD IDENTITY 



The Kubo formula for the conductivity describes the way external electromagnetic potential generates the motion 
of the charges. One needs special care in using this formalism when the charges in question are treated in the second 
quantized formalism. The point is that one should not use approximations which violate the particle number. In fact, 
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an uncontrolled number of particles created by an approximation renders the transport coefficients divergent when the 
IR limit is considered, whatever weak error is committed in the approximation. It is easy to control the total number 
of particles by global symmetries. But the complications start when the particles carry a conserved charge which can 
take positive as well negative values. In this case the conservation of the algebraic sum of the charge is not enough 
any more to protect the transport coefficient against IR divergences. In fact, imagine an approximation scheme which 
though preserves the total electric charge exactly but leads to the creation of an uncontrolled number of neutral pairs. 
The conductivity will obviously be IR divergent in such a scheme due to the polarization cloud generated by the 
erroneous approximation. Charged particles are usually handled by gauge theories and the approximations which 
leads to the creation of neutral pairs without violating the charge conservation are those which are invariant under 
global symmetry transformations but violate local gauge transformations. 

The best way to make sure that a symmetry is preserved is to check if the equation of motion in the given 
regularization and approximation scheme remains symmetrical. In the derivation of the equation of motion, (10), 
then one uses infinitesimal shift of the dynamical field variables which corresponds to a symmetry transformation. 
The action remains invariant and any contribution to the equation of motion must come from non-symmetrical terms, 
i.e. non-symmetrical parts of the dynamics are the only sources contributing to such a special combination of the 
equation of motion, called Ward identity 

In our case we generate the shift of the field variables by an infinitesimal gauge transformation, 



(CI) 



Since the electromagnetic fields appear as an external source only without dynamics they do not participate in this 
transformation. Nevertheless they appear in the discussion because the transformation (CI) is equivalent with 



where 



Y7± r j fx±fl ~ fx 



(C2) 



(C3) 



as long as the dynamics is gauge invariant. The transformation (C3) induces the change 



SS - 



-aV^V+e^-O^-^ + Ote 2 ), 



5^ 


Oo,x 


v> = 




O 

^],x 




5^ 


^J,X 


*!> = 



m 



-V+e x ^ ■O ,x-ip + O(e 2 ), 
-Vt Cx ^t • Oj, x ■ V + 0(e 2 ). 



(C4) 



Using this transformation as a reparametrization of the functional integral including the operators O m and with 
a corresponding source we find 



Ei ia T \ SWn \<j. a] \ „ 





a2 ~ \ 


[(» 


~h 2(Tj ' x ) 



5W [a,a] 1 8W [a,a] 

<J; 



3,x 



m 



'3,x 



8<Jj,x 



(C5) 



After integration in part, omitting the boundary terms and averaging over the disorder we have 

<SWo[o- - iv,(Tj,d-]~ 




ia T 



1 + — {<tq,x - i»x) 



S Wo [do — iv, <Jj , cr] 1 SWq [uq — iv, &j , cr] 



i>j,x rn ' 5dj tX 

By taking the functional derivative of the 'equation of motion' for the impurity averaging, 



(C6) 
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= / D[u]e 29 Jx - U / i — ■ J - — + -v x Wo[a +v,aj,a] 



we can write find the Ward identity 
= V r 



ia T \5W[a,a] ga T fS 2 W[a,a} 
1 + ~1\~ <T °' X 



1 a 2 _ \ <W[ct,ct] 1 <W[<7,<t] 

1 + TjO'i.a; I ; &i 



1pj,x 



m 



1p0,x h Jt V'O.t.xV'O, 

J 

By assuming the form 

W[a] = W + a-p*-^a-D-a + ^W^-a^a^ + 0(a 4 ) 
the first functional derivative of (C8) with respect to a taken at a = a = yields the conserved currents 



Jo,, = b« y % - + ^ [ W(o, tl x),(o,x),(«, y) , Jj,x = + ^s v , x p* 



(C7) 



(C8) 



(C9) 



(CIO) 



for arbitrary k and y. For U.V. finite model or when the U.V. divergences are logarithmic only the terms 0{a T ) can 
be neglected. 



APPENDIX D: FUNCTIONAL DERIVATIVES 

We give few details of the computation of the second functional derivative of the effective action (40) . The effective 
action is supposed to be the sum of terms like 



l[p\ = / V Ap (p x )d A p p . x . 

J X 



and the first four functional derivatives are the following 
Sj[p] 



5p a ,a 
Spa,a,8pp,b 
5pa,a5p0,bSpy,c 

5p a , a Sp0,bSp 7 ,c5pe,e 



I [S a , z d Pa V A , p d A p p , z + V A , a d A 5 a , z ] 

J z 

/ [Sa,zSb,zd p0 d Pa VA^d A p^ z + S a ,zd Pa VA,0d A 5b,z + 6b 1 zdp /3 VA, a d A 5a,z 
J z 



-I 



5a,z5 b ,z5c,zd Pl d Pl3 d Pa V A ^d A p^ z + S a!Z S b , z d P/3 d Pa V Atl d A d c< 



+Sa,zSc,zdp y d Pa V A ^d A 5 b , z + 6b, z 6 CiZ d Py d Pf) V Aia d A 5 a!2 



J. 



Sa,z8 b ,zSc,zS e ,zd Pe d Py d Pfl d Pa V A , p d A p p , z + 6 a , z 6 b , z 5 CtZ d Py d P/3 d Pa V Ate d A 6 etZ 
+Sa,z5 b ,z5 e ,zd Pc d Pfj d Pa V A ^d A 5 c , z + S a , z S c . z S e ^ z d P€ d Pl d Pa V At0 d A S b , z 
+Sb, z S c ,zS e ,zd Pc d Pl3 V A , a d A 5 a ,, 



We write V Ay(J = p K V KiAi(7 and find 



d Pa Va,o- = y a . A ,a + p K d Pa V KtA ,o- 
9 P0 d Pa V AtC7 = d Pfl V a . A ^ + d Pa Vp,A,o + p K d PI3 d Po y KtAta 
dp y dppdp«VA,<T = d Pi d Pfj V a . Ay(T + d Py d Pa Vp^ a + d P0 d Pa Vy t A,a + p K d Py d P0 d Pa V KiA ^ 
9 Pc dp y d P0 d Pa VA,a = d P€ d Py d P/3 V aiAi „ + d Pc dp^d Pa V^ At(r + d P€ d Pl3 d Pa V liAia 

+ d P 1 d P(i d P<> V e,A,a + Pudpc d p _ f d pf} d pa V K<A ,a 



(Dl) 



(D2) 



(D3) 



and 
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5p c 



l r 

-= [5a,zV a ,A^d A p (T ^ + S a!Z p K dp a V K! A,ad A p a ,z + P K V K . A , a d A 5 a ,z] 
a J z 



5 2 1[P\ 



: / \s a , z V a> A,ad A Pa,z + S a ,zp K d Pa V K , A , a d A Pa . z + {-l) lAl 5 a , z d A (p K V K> A, a ) 

Sa.zSb,z{d Pll V a ,A.<7 + d Pa Vf3, A .a + PkQ p.d Po V K . A ,cr)d A Pa ,z 
0Pa, a 0p/3,b J z _ 

+5 a ,z{V a , A ,f3 + p K d Pa V K ^f3)d A 5^ z + S b , z (Vp iA ,a + p K d P/3 V K ^ a )d A 5 a . z 
" ' - f 5a,z5b,z5 c , z {d Pi d Pfi V a .A,cj + d p ^d Pa Vp t A,<? + d Pl3 d Pa V 1 ,A,<j + p K d p ^d P/3 d Pa V Kt A,a)d A Pa,z 

J z . 

+(S a , z Sb,z{ d P0 V c,A,i + d p a v P,A,i + PKd P[i d Pa V K .A n )d A S c . z + perms.) 
/ Sa.zSb^Sc^Se^idptdp^dppVa^AM + d Pe d p _ t d Pa Vfj. A ,a + d Pc d p? d Pa V 7 ^ A ,a + d p ^d Pf) d Pa V e ,A, 

J z . 

+p K ,d Pt d p _ i d Pl3 d Pa V K ,.A,a)d A p a ,z + {5 a , z S b , z 5 CiZ (d Py d Pf} V a! A,e + d Py d Pa V ,A,e 
+d Pf3 d Pa V JiA ,e + p K d Py d Pll d Pa V K ,A,e)d A 5 e . z + perms. 



Sp a ,aSpf3,bSp 7 ,c 
Sp a ,a$Pf3,b$P~/,cSpe,e 



(D4) 



where \A\ = |^o| H • For homogeneous configurations p a-x = p a we find 



h[p] _ 

6 Pa, a 



= 



5 2 J[P] 

Sp a ,aSpf3A 

s 3 i[p] 

5p a , a 5pf3,b5p 7 , c 

s_Hp\ 

5p a ,aSpf3,bSp 7 ,cSpe,e 



or in Fourier space 



Sp a , P S 

sHp] 



Sp a ,p$Pf3,qSp~f,, 

5p a ,p5pp, q 5p 1 , r 5pe, ! 



S a ,z(V a .A,p + PKd Pa V K .A,p)d A 6 b ,z + Sb,z(Vf3,A,a + p K d P0 V KtA , a )d A 6 aiZ 

/ 5 a , z Sb,z (d Pl} V at A,-y + d Pa V p . An + p K d Pf1 d Pa V KyAn )d A 5 CyZ + perms. 
J z L 

= / S a!Z S biZ S c<z (d p ^ d Pli V ai A,e + d^d^Vp^e + d P i 3 d Pa V 1 ,A,e + p K d p ^d Pf} d Pa V K , A ,eW 

J z . 

+perms. 

= f e^+ lb H a . z 5b, z [(V a , A ,0 + P.d Pa V K , A ,p)(tq) A + (Vp, A ,a + p K d P ,V K , A , a )(ip) A ] 

J a.b.z 

= Sp+q.O [(Va,A,P + Pnd Pc ,V K ^(3){iq) A + {Vp,A,a + P K d pi3 V K ^a){ip) A ] 

= 5 P +q+rfl [{dppV^A^ + d Pa Vf). An + PKd Pfj d Pa V Kt A n )(ir) A + perms.] 
= 5p+ q +r+s,n {d Py d P0 V at A,e + d Py d Pa Vp,A,e + d P/3 d Pa V 7 ,A,e + PKd Py d Pll d Pa V K ,A,<d(,is) A 
+perms 



(D5) 



(D6) 



The first four functional derivatives of the full effective action on homogeneous background are given by 



S 3 T[p] 



'V. ; .,'V'- r ''''''' 



Sp a ,pSp/3, q 5p 7 , r 5p^ s 



-pct,/3,7,e r 
L p,q,r,s u p+g+r+s,0) 



(D7) 



25 



where 

= V a ,p(q) + p K d Poc V K ,p(q) + Vp, a (p) + p K d Pfi V K . a (p) ~ d Pa d pp U 
= V a .j3(q) + V{3. a {-q) + p K d Pa V Ki p(q) + p K d Pf3 V K , a (-q) - d Pa d P0 U 

r p;g,r = ( d P0 V a,-/{r) + d Pa V fjn (r) + p K d Pfj dp a V Kn {r)+ perms.) - d Pa d Pl3 d Py U, 

= p K d p/3 d Pa V Kt7 (r) + d Pi3 V a ^(r) + d Pc V fjri {r) + p K d Pa d Pl V K ^{q) + d Pi V a ^{q) + d Pa V Jt/3 (q) 
+P K d p/3 d p ^V K!a (p) + d Pi V fj . a {p) + d p/3 V 7!Ct (p) - d Pa d Pi3 d Py U, 
T p£'r,f = d P-, d P^A s ) + d Pi d Pa V p , t (s) + d P pd Pa V 7 , £ (s) + p K d Py d P0 d Pot V K , e (s) + perms. - d Pa d P0 d Py d Pe U 
= p K d Pl3 d Pa d Py V Kt e{s) + d Pa d Py Vp. t (s) + d Pfj d p _ t V a . t (s) + d Pl3 d Pa V lte (s) 
+ (Pnd P/1 d Pa d P€ V K ^(r) + d Pa d Pt Vp,~ ( (r) + d PlJ d Pt V an {r) + d P0 d Pa V en (r) 
+ {p^d Pfj d Pi d Pc V^ a {p) + d Pi d Pe V f 3, a {p) + d P0 d Pe V 7!CC {p) + d Pfi d Pi V t . a (p) 

+ (Pnd Pa d Pn ,d Pe V Kj 3(q) + d Pi d Pt V aJj {q) + d Pa d Pe V Jt p{q) + d Pa d Pl V t ^{q) - d Pa d p? d Pl d Pt U. (D8) 
and the dependence on p is suppresed. Therefore the second functional derivative up to 0(Sp 2 ) is 

Sp a , P Spp, q P+q '° p ' q p > q 2 l ™ +U[dp) 

J r 

T pT P = l Wr+-,0^ 7 ,r^,.r^y. (D9) 
J r,s 

The appplication of these formulae for the effective action (42) yields the second derivative 

r (o)a,/j = /r* + ±vSr*(p) + fv?r*(g) Jp ir-<( P ) + ^q±v;^ ^s (< z ) \ 

l v« - \ jq Ir st (< Z ) +J pivOr^(p) \T T {q)T q + \Y L {q)L q + \T T {p)T p + \T L {p)L p ) 
where T* = -d$U, T s = -d 2 U. The propagator is 



with 



°p+qfl L P ,q — I x y!jr ta i „ , r'+v?r" T ) 



q 2 



r D = (r* + v?r tt )r I - + ^-(v§r ts ) 2 . (D12) 



The third and fourth derivatives are of the form 

r°;°;° = v^ ,o(p) + v^Vb, (<z) + v£vb, (r) + a r* 

= ^r*(p) + \wlT tt {q) + l -V\T*{r) + d Q T* 

r°;°;f - v^ ,/(r) + $>v<,o(p) + d v ifi ( q ) + ^r* 

= ^vir ts (r)r, + ^ r s *(p) P , + ^ r st (<?)^ + a,r* 
r£;^ = *Vfc,/(r) + 9o^, fe (g) + <5 M 3 r s 

= ia r T (r)T r fe ^ + Ia r T ( g )T^ + id r L (r)z^ + h^T^L 1 * + s k , e d T* 
r°;°;°; s ° = v§^,o(p) + v 2 y ,o(<z) + v§vb,„(r) + v^.oW + a 2 r 4 

= ^v 2 r«(rt + iv 2 r%) + iv 2 r«(r) + iv 2 r»( s ) + a 2 r* 
r°;°;°; s fc = v 2 T/ , fe (s) + a 2 T4, (p) + d 2 v kfi ( q ) + d 2 v kfi (r) + a a fc r* 

= ^V 2 r ts ( S ) S fe + ^ 2 r st (p)/ + l -d 2 T st (q)q k + l -d 2 T st (r)r k + d^T* 

r°;° 9 ; r fc ;f = d 2 v k , e (s) + d 2 v e , k (r) + s k > e d 2 r s 

= \d 2 T T {s) T y + \d 2 T L {s)L k / + \d 2 T T {r) T y + \d 2 T L {r)L k / + 5 k '%T° . (D13) 
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APPENDIX E: ONE-LOOP W[a] 



This Appendix we briefly record the expressions needed for the computation of the Green functions occuring in 
the one-loop approximation of the generator funtional VF[cr]. This generator functional is non-trivial even for free 
electrons due to the exchange interaction. 



1. Densities 



We start with the expresions for the densities 



2 
and 



P* = - v G *,x,t+5,x + c.c. = I 9(-£ q ) (El) 



q 



Y^P*j = ~y X! G *' x > t+l5 > x +J + c - c - = J B(-£; q )cosgia (E2) 



1 ih 

p = 



3 

3 J>0 

in terms of the electron propagator G p . The local potential is given by 

n 



Vi 



TrlogL*- 1 = -i / log[£-£; q -He q ]. (E3) 



Charge conjugation acts as p — > p c = 6h 2 /ma 2 — p, p* — ► p c * = B p — p* and p* — > p c * = p* where the band 
contribution is 



/ 1 (E4) 

Jq 



2. Two-point functions 



We need 



£0'° = -in 



J G r Gp-\- r , 

D / = -— e" <2 £ / G r G p+rS ma( rj + 
' ma J r V 2 / 



Df = 



z^ 3 



^2 ^ G r G p+r sin a (rj + y- ) sin a (r k + y ) , (E5) 



where the Fourier transform is defined as 



^ = / ^e*^. (E6) 

It is straightforward to perform the energy integrals. The results are shown for the external momentum p = (0, 0,p) 
after having carried out the replacement p — ► p — <7q , 

D°>°(oj,p 2 )=2V f ® ( ~ E i AE T 2 +i* I e(-£? q )0(£; q+ p)[*(A£; q ,p-M + *(AS q ,p + M] 

D°J( W ,p) = -^P / e (-^) sina (g3 + f) 
am 7 q (A£ q , p ) 2 - ?i 2 w 2 

y e(-S q )e(^ q+p ) sin a ( 0j + y) [<5(A£ q , p - Rw) - <S(A£ q , p + Ml 



ihir 
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L> 0j (-w,-p) = -D^{-u,p) = D j >°(uj,p) 



+ / e(-£ q )9(.E q+p ) sin 2 a gi [<5(A£; q , p - ftw) + <S(A£ q , p + ftw)] 

a m ./q 

2ft 2 /■ 0(-£? q ) sin 2 a (qs + §) A£ q , p 



£>,p 2 ) 



(A£ q , p ) 2 - ft 2 co 2 

--^ / e(-£ q )e(£ q+p )sin 2 a ( gi + I) [<5(A£ q , p - ftu,) + 5(AE^ p + ftu;)], (E7) 



where A_E qiP = _E q — i? qp and V denotes the principal value prescription. The singularity and the Dirac delta in 
the integrand render these expressions unpractical for numerical evaluation. We used Wick rotated integrals and 
considered the dependence in Euclidean energies in the numerical work. The corresponding expressions are 



7 q (A£ q , p ) 2 + ft V 



£°'V,P 2 ) 



~ QJ 2iuh 2 f Q(-£ q )sina( g3 + f)l 

L> ' J (bJ,p) — / s 

am 7 q (A£ q , p ) 2 + ftV 



■ rp 2 2ft 2 /" 9(-£; q )sin 2 agiA J B, 



^ T (^,P 2 ) 



2ft 2 y 

a 2 m 2 J q 



(A£ qjP ) 2 + ft V 
e(-£ q )sin 2 a( g3 + f) A£ q>p 



These expressions show clearly the need of keeping a non-vanishing energy ftw in any computation carried out in 
a finite system. In fact, the limit u> —* generates a Dirac delta in the integrands which yields singular results in a 
finite system where the density of states is the sum of Dirac delta peaks. 
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